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Abstract. We discuss recent evidence that currently accepted mass- loss rates 
may need to be revised downwards, as a consequence of previously neglected 
"clumping" of the wind. New results on the radial stratification of the cor- 
responding clumping factors are summarized. We investigate the influence of 
clumping on the ionization equilibrium of phosphorus, which is of major rele- 
vance when deriving constraints on the clumping factors from an analysis of the 
FUV Pv resonance lines. 



1. Introduction 

While our understanding of the outflows from luminous OB-stars was thought to 
be well established, recent evidence indicates that currently accepted mass- loss 
rates may need to be revised downwards by as much as a factor of ten. This 
is a consequence of previously neglected "clumping" of the wind, which affects 
mostly those diagnostics which are sensitive to the square of the density, p (such 
as recombination lines or free- free continua). 

Considering that numerous stellar-evolution calculations have demonstrated 
that changing the mass-loss rates of mas sive stars by even a factor of two has 
a dramatic effect on their evolution (e.g., Mevnet et alJll994h . it is evident that 
such revisions would have enormous implications, not only regarding evolution, 
but also regarding the feed-back from massive stars. In this article, we will 
summarize the knowledge which has been accumulated lately and consider the 
question concerning the REAL mass-loss rates from massive star winds. 



2. Standard mass loss diagnostics of luminous OB-stars 

Traditionally, the mass-loss rates, M, of luminous OB stars have been inferred 
from, primarily, three types of measurement (see also de Koter, this volume): 
1. The strengths of UV P-Cygni profiles were the first diagn ostics to be used to 
measure wind-densities. This approach has been pioneered bv lLamers fe Morton! 
(1976) in their work on the "Mass ejection from the 04f Star Zeta Puppis", 
where they derived a mass-loss rate of 7.2 ± 3.2 • 10 -6 MQyr -1 , a number which 
is still valid! The most c ommonly used met hod to analyze P-Cygni lines is 
the so-called SEI- method ( Lamers et al.llT987 l. which has been firstly applied 
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bv iGroenewegen Sz Lamersl (Q.989). Though most of the UV P-Cygni lines are 
resonance lines from dominant ions being linearly dependent on density and 
thus remaining uncont animated from direct clumping effects (but see Sect. 17^1 . 
only the product Mq can be inferred from these lines, where q is the ionization 
fraction of the corresponding ion. Moreover, most of these lines are inevitably 
saturated in stars with strong winds, since they arise from abundant elements, 
and only lower limits on M can be provided in these cases. 

2. Thermal radio and FIR continuum emission samples the outermost and in- 
termediate region of the wind (0(100, 10 R+)), r e spect i vely. The corresponding 
methods ba se on the work by IWright fc Barlowl (|l975h / IPanagia fc Fellil (| 19751 ) 
(radio) and lLamers fc Waters! (|l984allbD fIR- excess) . They have been applied, 
e.g., bv lAbbott et al.l (|l98lf )~ and|L amers Leithererl (|l993l ) to derive OB-star 
radio mass-loss rates, and bv lLamers et al.l ~ 1984T ) to investigate the IR excess 
of £ Pup by means of IRAS observations. Since the involved processes are dom- 
inated by free- free and bound- free transitions which scale with p 2 , inferences 
from the measurements are extremely sensitive to clumping in the wind. 

3. H Q emission, usually modelled using NLTE atmosphere codes, also depends 
on p 2 , and is therefore also sensitive to clumping, but samples the innermost 
portion of the wind ( < 2R+). The ide a to ex ploit H Q as a standard mass-loss 
indicator goes b a ck to Klein Castorl (ll97ST). and has been firstly applied by 
iLeithererl (|l988l ). iDrewl Il990) and lLamers Leitherer l (^£93). Further refine- 
ments have been provided bv lPuls et al.l (|l996h . 



3. Clumping — some basic facts 

The present hypothesis states that clumping is a matter of small-scale density 
inhomogeneities in the wind 1 , which redistribute the matter into clumps of en- 
hanced density embedded in a rarefied, almost void medium. The amount of 
clumping is conveniently quantified by the so-called clumping factor, f c \ > 1, 
which is a measure of the over-density inside the clumps (compared to a smooth 
flow of identical average mass- loss rate). 2 As already pointed out, diagnos- 
tics that are linearly dependent on the density are insensitive to clumping, 
whilst those sensitive to p 2 will tend to overestimate the mass- loss rate of a 
clump e d wind, by a fact o r \fJ7\ - For further de tails, see, e.g 
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ciumpeu w mu, uy a raci o r y ; r i . ror runner ueians, see, e.g., auooll 
( 19811 ). lLamers fe Waters! (Il984bl ). ISchmutd (|l995f ) and lPuls et all 12006) 



Until to date, the most plausible physical process responsible for small-scale 
structure formation in massive star winds is the so-called line-driven instability, 
found already in the first time-dependent hydrodynamical simulations of such 
winds (lOwocki et alJ Il988f ) ; for recent investig ations, see iRunacres Qwockil 
1 20021 . 12005ft with respect to 1-D simulations and lDessart fc Qwockil (|2003l . 12005) 



for first attempts to include 2-D processes. 

Though predicted from the earliest hydro-models on (and even before, see 
iLucv &z Solomonlll970l ). it took some while to incorporate clumping into the at- 



1 in contrast to large-sca le inhomogeneities such as co-rotating interaction zones (e.g., 



in contrast to large-scale inhomogeneities si 
ICranmer fc Ow ocki 1996' and references therein) 



2 An alternative description bases on the volume-filling factor, /v = f c] 
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mospheric models of massive stars, firstly for Wolf-Rayet (WR) atmospheres, in 
order to explain 

(a) the strength of the observed electron scattering wings (/>depe ndent) in par - 
allel with the strength of the underlying (/9 2 -dep.) emissi on lines (lHillierlll99]] 'l. 

(b) the so-called momentum problem of WR winds (e.g.. ISchmutz 19951 ). 

(c) the presence and variability of sub-structures in WR emission lines, by in- 
yoking supersonic turbu lence leading to clumping factors of the order of f c \ ~ 9 
(Moffat fc Eobertlll99l N l. 



4. No clumping in OB-star winds? 

In contrast to the case of WR winds, the diagnostics of OB-star winds did not 
render any necessity for (significant) clumping until recently. In particular, two 
major arguments supported the view of a rather smooth, almost unclumped flow: 



1. The investigations bv lLamers Leithererl ( 19931 ) and lPuls et al.l (|l996h showed 



that Hfy and radi o mass-loss rates are similar for a large sample of stars (but see 
also lDrewl (|l990) who noted a discrepancy by a factor of two, the former being 



larger). Since H a forms in the lower and the radio emission in the outer wind, 
this would imply a similar degree of clumping in both regions, which seemed to 
be improbable and is also in contradiction to theoretical predictions. 
2. The observed wind-momentum rates were found to be in rather good agree- 
ment with theoretical pred ictions for this quantity, obtained from different , 
independent investigat ions flVink et al J 1200(1 . iKudritzkil 120021 . IPuls et al J 120031 
iKrticka fc KubatJ I2004D . A pure coincindence of this agreement seemed to be 
rather unlikely. 

Taken together, it was concluded that clumping effects in OB-star winds should 
be negligible. 



5. Indications of (significant) clumping in OB-star winds 

As stated in the introduction, recent evidence points to the notion that this 
conclusion might be incorrect. In the following, we will summarize the different 
indications for significant clumping in OB-star winds accumulated so far. From 
the observational side, there is, to our knowledge, only one direct evidence for 
clumping. From a temporal analysis of Hell 4686, lEversberg et al. | (119981 . m 



particular their Fig. 2) found "outward moving inhomogeneities" in the wind of 
C, Pup, from regions near the photosphere out to 2 R+. 

All other evidence is indirect, and the derived clumping factors cover a large 
range. 

(a) From polarimetry and using simplified models, clumping factors of the order 
of /d=2 are found (Davies, this volume). 



(b) The an alysis of radio and submm data from e Ori and £ Pup ( Blo mme et alJ 



12001 12003) indicates a clear submm excess in both cases, which can be explained 
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by (enhanced) clumping in the intermediate wind (~ lO-R*) (see also Blomme, 
this volume). 

(c) NLTE model atmosphere analyses of UV spectra (partly incl. the optical 
range) of various O-stars indicate clumping factors of the order of f c \ = 10. . .50, 
but show that only few lines are suited to discriminate between clumped and 
unclumped flows. From these analyses (cf. Table P), it was concluded that 
clumping, if present, should begin rather close to the wind base, at a few tens 
of kmsec^ 1 . This finding is in contrast to hydro dynamical simulations, which 
show that the line-driven instability needs a certain time to grow and to become 
non-linear, so that significant clumping is expected not before 1.2 to 1.3 R±. 



Table 1. Evidence of clumping from UV diagnostics (NLTE). 



authors 


objects 


indicator 


fcl 


comments 


Crowther et al. 


AV232(07Iaf+) 


Pv 


10 


other lines barely affected 


(2002) 


SMC 






by clumping 


Hillier et al. 


AV83(07Iaf+) 


Pv/strong 


10 


if clumping is important, 


(2003) 


SMC 


UV photo- 




it must begin at relatively 






sph. lines 




low velocities (30 km/s!) 


Bouret et al. 


SMC dwarfs 


Ov 


signi- 




(2003) 






ficant 




Bouret et al. 


HD 190429A (04If) 


Pv, Ov, Niv 


25 


reduction of M by factors 


(2005) 


HD 96715 (04V((f)) 




50 


of 5 and 7 










clumping must start 










at the wind base 



(d) From d etailed investigati o ns of the wind-momen tum rates of a large s ample 
of O-stars, iPuls et alJ (|2003h . iMarkova et aD (I2004T ) and lEepolust et al.l (120041 ) 



found that supergiants with H a in emission lie above the theoretical wind- 
momentum luminosity relation (WLR), whereas the rest fits almost perfectly . 
Since the WLR should be independent of luminosity class (e.g.. lPuls et al.l ll996b 
this discrepancy was interpreted in terms of clumpy winds, with f c \ ~ 5, and 
mass-loss rates reduced by factors between 2 and 3. 



(e) A compelling, independent indication of clumping com es from SEI analy- 
ses of the P-Cygni P VAA1118, 1128 resonance line doublet ([Massa et al . 2003, 
iFullert on et al.l 120061 . see also Massa, this volume), which has only become 



widely accessible since the launch of FUSE. Because phosphorus has a low cos- 
mic abundance, this doublet never saturates in normal OB stars, providing useful 
estimates of not only Mq (cf. Sect. but also M itself, when P 4+ is the dom- 
inant ion - as it is implied at least for mid-0 star winds. These mass-loss rates 
turned out to lie considerably below those inferred from other diagnostics such 
as H Q or radio emission. The most reasonable way to reconcile these results is to 
invoke extreme clumping in the wind (Pv as a resonance line remains insensitive 
to direct clumping effects), indicating clumping factors of the order of f c \ « 100 
(or even more). Accordingly, the actual mass-loss rates should be much lower 
than previously thought, by a factor of ^> 10. Further comments are given in 
Sect. El 
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(f) Using similar methods. iPrinia et al.l ( 2005h showed that the unsaturated P 
Cygni lines in lower luminosity B supergiants give, again, a f a ctor o f 10 lower 
mass-loss rates than theoretically expected. ICrowther et al. see also 

Crowter, this volume), on the other side, found reasonable agreement between 
observed and predicted mass- loss rates for early /mid Bla supergiants, with simi- 
lar (though unconstrained) clumping factors in the lower and intermediate wind. 



6. A combined Ha/IR/mm/radio analysis 

Whereas most of the above investigations are concerned with a global (i.e., ra- 
dially almost constant) clumping factor or derive this quantity for certain wind 
regions only, there is, of course, the additional qu estion regarding the radial 
stratification of f c \. To this end, iPuls et all ()2006l ) recently performed a self- 
consistent analysis of H a , IR, mm and radio fluxes, thus sampling the lower, 
intermediate and outer wind in parallel. This study comprises a sample of 19 
Galactic O-tvpe supergiants and giants wit h we ll known stellar param eters (from 
iMarkova etahl |2004 iRepolust et aDl2004l and lMokiem et al.ll2005h . employing 
own measurements/archival data for H Q , IR/mm fluxes (SCUBA) and new VLA 
observations, including objects with H a in absorption, i.e., low density winds. 

A major result of this investigation is that in weaker winds the clumping 
factor is the same in the inner (r < 2R+) and outermost regions. However, 
for stronger winds, the clumping factor in the inner wind is larger than in the 
outer one, by factors of 3 to 6. This finding points to a physical difference in 
the clumping properties of weaker and stronger winds (and may be related to 
the excitation mechanism of the structure formation). In terms of mass-loss 
rates then, we find M(radio) ~ M(H a ) for weak winds with H a in absorption, 
whereas for all stars with H a in emission we obtain M(radio) ~ 0.4. . . 0.5 M(H a ), 
consistent with th e arguments given by Markova et al./Repolust et al. and the 
earlier findings bv lDrewl (|l99(f l. 

Unfortunately, this analysis is hampered by one severe restriction. Since all 
diagnostics employed have a p 2 dependence, only relative clumping factors could 
be derived, normalized to the values in the outermost, radio-emitting region. In 
other words, M(real) < M(radio), since until now the clumping in the radio 
emitting region is still unknown. Only if f c \ (radio) were unity, we would have 
M(real) = M (radio). Thus, the issue of absolute values for M still remains 
unresolved. 



7. The PV problem 

The major result of the investigation bv lFullerton et all (|2006T ) (see Sect. Op) is 
displayed in Fig. ^ In this figure, they have plotted, as a function of T e fj and 
luminosity class, the quantity 

q t = (Mq) ohs ^ (Mg)obs ^ _(?)_ (1) 

^Ha/radio My/fe vTd' 

for P 4+ , where angle brackets denote spatial averages. The quantity in the nom- 
inator has been measured from the Pv lines, and the corresponding H a /radio 
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T e „ [kK] 



Figure 1 . Derived estimates for c w, (Eq. U) of P^+, as a function of T e g 
and luminosity class (for details, see iFullerton et aI1 l2006). Overplotted are 
predicted values for this quantity, as obtained from own calculations. Bold: 
unclumped models; other curves: clumped models, with / c i= 9, 36, 144 (dot- 
ted, dashed, dashed-dotted), respectively. 



mass-loss rates have been taken from the literature. If the winds were un- 
clumped, (/est would correspond to the average ionization fraction of P 4+ , whereas 
for clumped winds this quantity is modified by f^ ' 5 - Fullerton et al. now argue 
that P 4+ should be a dominant ion at mid O-type, and in this case the derived 
clumping factor would be f c \ = 0(100) at T e fj« 40,000 K. Additionally, how- 
ever, they report on test calculations performed with unclumped models which 
show that, on the contrary, P 4+ should become a dominant ion only below 07. 
If this were true, Fig. ^would imply f c \ ~ 10,000 in this temperature regime! 

Since this is VERY unlikely, we have investigated the influence of clumping 
on the ionization fraction of P 4+ . Due to the enhanced density inside the clumps, 
stronger recombination is expected (this is the indirect ^ effect of clumping), which 
should change the picture (see also iBouret et 8,1.1 12005). By means of our model 
atmosphere code F ASTWI ND (Puis et al. 2005) and using the phosphorus model 
atom provided by iPauldrach et alJ ~ 200 lT ). we have calculated a sequence of 
clumped models with M^Tci = const, i.e., models which have identical p 2 - 
diagnostics mass-loss rates. As shown in Fig. [51 increased clumping indeed 
shifts P 4+ as a dominant ion towards higher T e fi: for unclumped models, it is 
dominant at 08/7, whereas for the largest clumping factors it dominates at 05. 

Using these models then, the "observed" run of g e st can be reproduced with 
highly clumped models (/ c i= 144, see Fig. for almost all luminosity classes 
and except for the very hottest temperatures (see below). Thus, our calculations 
confirm the hypothesis stated by Fullerton et al., i.e., clumping is possible to 
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Figure 2. Ionization fractions of phosphorus, as a function of velocity and 
spectral type, for supergiants. Black: Unclumped models (bold - Pv, dotted 
- PlV, dashed - Pvi). Grey: Ionisation fractions of Pv, for clumped models 
with fd— 9, 36, 144 (bold, dotted, dashed), respectively. 



explain the observations, and mass-loss rates might indeed be lower by factors 
of 10 or even more. 



8. Implications 

As we have seen, there seems to be a (physical) difference between thinner and 
thicker winds. For thinner winds, there is a similar degree of clumping in the 
lower and outer wind, whereas for thicker winds clumping is stronger in the 
lower part. The actual mass-loss rates depend on the clumping in outer wind, 
which is still an unsolved issue. 

If the outer winds were unclumped, our results would be consistent with the- 
oretical WLRs. In this case then, one would meet a severe dilemma with the 
results from the F(UV), which might hopefully be expla ined by additional ef- 
fects from X-rays emi tted due to clump-clump collisions (Fel dmeier et al.lll997l . 



iPauldrach et al.ll200ll ) . X-rays might also help to solve the problem encountered 
for Pv at highest T e g. 

If on the other hand, the (F)UV values were correct, the outer wind must be 
significantly clumped, and the present match of "observed" and predicted WLR 
would indeed by only coincidental. This scenario would imply a number of severe 
problems, not only for radiation driven wind theory, but, most importantly, 
concerning the stellar evolution in the upper HRD and related topics. A possible 
way out of the latter problem has been suggested by N. Smith (this volume). 
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